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Scanning tunnelling spectroscopy (STS) measurements find that the surface of Bi-2212 is charac- 
terized by nanoscale sized regions, "gap patches," which have different magnitudes for the (i-wave 
energy gap Ao(r). Recent studies have shown that the tunnelling conductance can be fit using a 
BCS-type density of states for a d-wave superconductor with a local quasiparticle scattering rate. 
The fit is made with a scattering rate which varies linearly with energy and has a slope Q;(r) that 
is positively correlated with the local value of the gap. We consider first the question of what is 
actually measured in such an experiment. To this end we revisit a model of quasiparticle scattering 
by impurities and spin fluctuations which was previously used to describe the lifetimes of nodal 
quasiparticles measured by angle- resolved photoemission (ARPES). We argue that the broadening 
of the local density of states is determined, except in the case of localized impurity bound states, 
by the imaginary part of the self-energy of the system averaged over a small region. The size of this 
region is set by a mean free path which depends upon the energy. At low energies, this region is 
found to be significantly larger than a gap "patch" , so that the density of states measured by STS 
is homogeneous in this energy range. At higher energies where the mean free path is comparable 
with the patch size, the density of states is inhomogeneous. We show that a local self-energy in 
the impurity-plus-spin fiuctuation model, while not strictly linear, yields a local density of states 
(LDOS) nearly identical to the full theory, and argue that it is consistent with the STS data as 
well as the phenomenological linear scattering rate extracted from experiment. We also explore the 
qualitative consequences of this phenomenology for the spectral widths observed in ARPES and 
predict the existence of Fermi arcs in the superconducting state. 

PACS numbers: 74.25.Jb, 74.20.Fg, 74.72.-h 



I. INTRODUCTION 



For many years, it has been recognized that the prop- 
erties of the quasiparticle states near the nodes of the 
(i-wave superconducting gap in the cuprates were quite 
different from those near the gap maxima or antinodes. 
Transport estimates of nodal mean free paths at low tem- 
peratures range from tens of jim in the clean stoichio- 
metric YBCO compoundsJl| to hundreds of A in dirt- 
ier materials like BSCC0[2]. ARPES measurements on 
BSCCO, with much cruder energy resolution, have re- 
ported residual nodal widths of the momentum distribu- 
tion curve (MDC) of about 10meV[3], which with Fermi 
velocity values of 10^ cm/s imply much shorter mean 
free paths, of order 30 A. While this discrepancy is not 
currently understood, both transport and ARPES nodal 
mean free paths are significantly longer than those ex- 
tracted from ARPES for antinodal quasiparticles. These 
scattering rates are of order 25-40meV at low T in the 
optimally doped superconducting state [1], correspond- 
ing naively to mean free paths of order a unit cell size. 
One explanation for this "dichotomy" between nodal and 
antinodal states that has been proposed is the existence 
of intense magnetic scattering near the antinodal point 
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in the underdoped materials due to the proximity of the 
wave vector Q connecting the antinodal points by an an- 
tiferromagnetic nesting vector. This scattering has also 
recently been related by several authors to the properties 
of the "arcs" of coherent Fermi surface in the pseudogap 
state [i, El, 0. 

Recently, a new perspective has been brought to bear 
on the properties of these quasiparticle states by STS. In 
Ref. [8], Alldredge and collaborators succeeded in fitting 
STS data on Bi-2212 at various dopings to an extremely 
simple form, that of the local density of states (LDOS) 
of a d-wave superconductor with a single position and 
frequency dependent broadening function or scattering 
rate F(r, cc;). 



uo -f ir(r, uo) 



(u; + ir(r,u;))2-A2(r)0d(k)2 



(1) 



Here 



■i- Electronic address: djs @physics.ucsb.edu| 



r(r,a>) = ri(r) +r2(r,cu) = ri(r) + a(r)a>, (2) 

with (/)(i(k) = {coskx — cos ky)/ 2 and a{r) a constant 
of proportionality which was found to vary from point to 
point over the sample surface, and to be strongly and pos- 
itively correlated with the size of the local gap parameter 
Ai(r). The frequency independent term Fi simulates the 
effect of near unitarity limit in-plane scatterers, which are 
observed as "native defect" resonances in conductance 
maps of the Bi-2212 surface. Values of F extracted from 
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this fit procedure tended to be much smaher than quasi- 
particle widths determined by ARPES, ranging at en- 
ergy uo = (Ai) from about 2 meV for an optimahy doped 
sample to about 20 meV for a strongly underdoped sam- 
ple. These results, and the apparent contradiction with 
ARPES, lead us to consider the following questions: 

• What sort of scattering rate is actually being mea- 
sured in the STS measurement? 

• When is it possible to think of a local probe mea- 
suring a scattering rate averaged over a region or 
over disorder, and what length scales define such a 
region? 

• How can the spatially modulated r2 at low energies 
implied by Eq. ([2]) be reconciled with the observed 
homogeneity of the low-energy quasiparticle states 
seen in STS? 

• Are the results of Ref. [8|] possibly consistent with 
scattering rates other than pure linear in cj, deriv- 
able from microscopic models? 

• How are the small values of the STS-determined 
scattering rates to be reconciled with ARPES? 

In the following we attempt to answer these ques- 
tions within a phenomenological framework which as- 
sumes as its starting point that there are nanoscale en- 
ergy gap inhomogeneities. Various suggestions regarding 
the origin of this "gap patch" inhomogeneity have been 
proposedl, [m, E E, Q [H. Here we will sim- 
ply assume that for T ^ Tc, the superconducting gap 
has such nanoscale patch inhomogeneities. We will then 
proceed in several stages. In Section |TT1 we discuss the 
local density of states of a disordered system, and the 
conditions under which a system may be considered to 
be locally self-averaging. We then extend these consid- 
erations to inelastic scattering, and discuss what may 
be deduced from the qualitative aspects of the STS re- 
sults. In Section [TTTl we introduce a model for scattering 
from impurities, including native defects in the Cu02 
planes, and dopant disorder away from the planes. To 
the disorder- averaged impurity self-energy we then add 
the self-energy due to scattering from spin fluctuations, 
and discuss the anisotropy of the total scattering in mo- 
mentum space. We then calculate the LDOS from the 
local Green's function in Section \IV\ and compare the 
results with those computed from a model "local self- 
energy". We show that the spectrum may equally well 
be fit by an assumed linear scattering rate, as proposed 
in Ref. [8i] . Finally, we discuss the extension of the results 
obtained for the homogeneous d-wave superconductor to 
the inhomogeneous case. We discuss how the spatially 
homogenous low-energy conductance seen in the STS ex- 
periments for CO < Aq and the heterogeneous behavior for 
00 > Aq may be understood within this framework. In 
addition, we point out the implications of Ref. [8] for the 
interpretation of ARPES spectra in the superconducting 
state. 



II. GENERAL CONSIDERATIONS 

We begin by assuming that STS probes the LDOS, i.e. 
the imaginary part of the local Green's function G(r, r); 
that is, we assume that the tunnelling matrix elements 
are constant. For the moment we focus on disorder in the 
normal metallic state, and neglect interactions. In this 
case one should in principle calculate G for a given config- 
uration of impurities. There is no notion of a spatially or 
disorder averaged self-energy which need appear, if one 
can calculate G(r, r) for a system with disorder. Indeed, 
the exact local Green's function may be written 

Oi.^:^^-i:J^. (3) 

n 

where V^n(r) are the exact eigenfunctions of the many- 
electron system in the presence of a given random disor- 
der configuration. Such a calculation would include all 
interference processes from these many impurities. 

Note that even if we are interested exclusively in the 
local Green's function, it is far from obvious that one 
may speak of a local self-energy or scattering rate, as 
proposed by the STS experiments. This is because the 
self-energy is formally a nonlocal quantity as is seen from 
the Dyson equation for G(r, r): 

G(r, r) = Go(0) + Go(r - r')I](r', r")G(r", r), (4) 

where Gq is the Green's function of the homogeneous 
system. S may be treated as local, I](r, r') ~ I](r)J(r — 
r') only if E(r, r') decays sufficiently rapidly. 

Under some circumstances is it appropriate to approx- 
imate the many-impurity Green's function by a homoge- 
neous disorder-averaged Green's function G. In a dirty 
metal, the wave functions tp{r) decay on a length scale of 
the mean free path £ of the system, as do the disorder av- 
eraged Green's functions and self-energies. This will then 
be a good description of the system provided there are no 
bound state wave functions present where electrons are 
confined to within a radius substantially smaller than 
i. Near the bound state energy (5(r, r; co) will provide a 
poor description of the spatial and spectral properties of 
this state. 

Similar considerations apply to quasiparticles in super- 
conductors. In general, a local measurement by STS will 
probe wavefunctions spread out over a mean free path. 
However, one should be wary in a highly anisotropic sys- 
tem like a d-wave superconductor of assuming the rele- 
vance of mean free paths extracted from other measure- 
ments. For example, transport measurements at low T 
yield the mean free paths of nodal quasiparticles only, 
and can be very long, whereas an STS experiment prob- 
ing the system locally at bias co will involve quasiparticles 
with a momentum spread uo/vp and these may well have 
shorter mean free paths. We will see that results of the 
STS measurements may be used to determine the typi- 
cal mean free path of the excitation being measured. In 
general, however, the system may indeed be taken to be 



3 



self- averaging and a treatment in terms of a disorder- 
averaged G should be adequate. 

An exception occurs again in the case of impurity 
bound states, known to occur in the cuprate materials for 
Cu substituents in the Cu02 plane. As observed by STM, 
in these cases a very narrow spectral resonance is found 
at frequency VLq on the impurity site and near neighbors. 
The naive calculation of the size of these states at reso- 
nance yields 



4 = 



^0 



02 



(5) 



where Aq is the gap maximum and (fo = '^F/(7rAo), ex- 
cept in the nodal directions, where ~ 1/r and di- 
verges. Note that the electron unbinds when the reso- 
nance moves into the continuum VLq > Aq, corresponding 
to weak impurity potential. More detailed analysis shows 
that the length scale ^5 given above represents the asymp- 
totic long-range decay of the wavefunctions in non-nodal 
directions; most of the decay away from the impurity site 
actually takes place over an atomic distance. In any case, 
for in-plane impurities, <C /, so they should be (and 
are) imaged individually by STS. On the other hand, the 
out-of-plane dopant atoms in the BSCCO-2212 material 
have weak potentials and no resonance, corresponding to 
an essentially infinite ^5, i.e. a disorder averaged calcu- 
lation of the local Green's function is again appropriate. 
We conclude that a disorder- averaged (5 is a good rep- 
resentation of the local exact G(r, r) over a length scale 
of size i except within a distance of strong impurities. 
Therefore the notion of a local self-energy or scattering 
rate will also have meaning, except for a small set of ex- 
ceptional points. 



III. MODEL FOR SCATTERING 

The purpose of this section is to present the results of 
a scattering-rate analysis of spin fluctuations and impu- 
rities which will be used to determine the local G. We 
assume that a self- aver aged description is appropriate, 
and investigate, in the BSCCO system, what the sources 
of scattering are, and how they enter the local G. We can 
then ask if a description in terms of a local self-energy or 
scattering rate is sensible or not, and whether or not a fit 
to the LDOS obtained from G can also be obtained with 
a linear scattering rate as proposed in Ref. Q . This will 
give us information as to the robustness of such a fit. 



A. Elastic potential scattering 

The elastic scattering in BSCCO can arise from im- 
purities and disorder inside the conducting copper ox- 
ide planes (in-plane impurities) or due to impurities in 
the neighboring metal oxide planes (out-of-plane impuri- 
ties) . While the in-plane impurities are well characterized 
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FIG. 1: (Color online) The first two terms of Eq. ^ (ro 
and the t\ contributions) to the elastic scattering rate as a 
function of kx and ky in the first Brillouin zone for different 
quasiparticle energies uj. From left to right and top to bottom 
ijj is chosen to be cj = O.lAo, = 0.2Ao, a; = 0.5Ao and 
uj = Aq. The red line shows the Fermi surface obtained from 
the single t — t' band described in the text. The inverse range 
of the impurity potential is n — 0.2, corresponding to strong 
forward scattering. The color scale is the same for all four 
plots. 



by a point-like isotropic scattering potential, the out-of- 
plane impurities are better described by a long range e.g. 
Yukawa type impurity potential leading to a more or less 
pronounced forward scattering [l^, Since the 

out-of-plane impurity scattering is relatively weak it can 
be treated within the self-consistent Born approximation. 
The elastic self-energy is then given by 

S,i(k,c^) =n,^|F(k,kO|V3G(k^c^)r3 (6) 

k' 

where Tq, are the Pauli matrices in particle hole space, Ui 
is the impurity concentration and V(k, k') characterizes 
the momentum dependent impurity potential. Within 
our calculations we will model it by a screened exponen- 
tial falloff of the form 



|y(k,k')| 



/|2 



(7) 



where characterizes the range of the impurity po- 
tential. For n ^ 00 the case of isotropic scattering is 
recovered. The full Green's function G(k, uo) can be cal- 
culated from the Dyson equation 



G-^(k,t^)=Go-nk,c^)-S,i(k,t^) 



(8) 



where Go(k, a;) describes the pure superconducting state 
without scattering. Decomposing the self energy into its 
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different Nambu components leads us to the following 
three equations which must be solved self-consistently 

Sei,o(k,a;) = ^ |V-(k,kO|^ ~2_~2 _X2 ^ 

Sei,i(k,u;) = -n.^|y(k,kO|^ ,^ _ (10) 
and 

Sei,3(k,c^)=n.^ink,kO| ^^ ~, (11) 

- 4' - ^k' 

Here = - Eei,o(k, cj), Ck = ek + Sei,3(k,ct;) and 
= Ak + Sei,i(k, ijj) are the renormalized quasiparticle 
energy, band structure and pairing potential respectively. 
For our numerical calculations, we will use a simple tight- 
binding parameterization of the band structure including 
nearest and next nearest neighbor hopping 

£k = — 2t (cos krj. + cos ky) — ^t' cos cos ky — \i (12) 

with a d-wave order parameter 

Ak = {coskx — cos ky) . (13) 

To model the Fermi surface in the cuprates we use 
= — 0.35t and fi = — l.lt, with energy measured in 
units of nearest neighbor hopping t. The temperature 
dependence of the gap is parameterized by the simple 
form 



Ao(T) = Aotanh J^a^ ^ - ij (14) 

with a = 3, 2Ao/kTc 6, and Aq 0.2t. 

The elastic scattering rate, broadening the quasipar- 
ticle state with energy a;, is determined by the denomi- 
nator of the full Green's function and can be calculated 
approximately from the pole of G as 

rei(k,cj) = -Im ^Eei,o(k,cj) + ^i;ei,i(k,cj) 

+ 5iSei,3(k,C^)) (15) 
CO / 

For quasiparticles at the Fermi surface only the imag- 
inary parts of Eei,o and Sei,! contribute to the elastic 
scattering rate. To compare our results for different im- 
purity ranges, we have normalized the scattering rates in 
a way that the (normal conducting) scattering rate for Tc 
taken at the nodal point of the Fermi surface Icat equals 
O.IAq. All calculations for the elastic as well as for the 
inelastic scattering rate have been performed for a low 
temperature T = O.lTc unless otherwise stated. 

In Fig. [T] we show the imaginary part of the tq and 
the Ti contribution of the elastic self energy to Fqi = 




FIG. 2: (Color online) Same as Fig. (TJ but with less pro- 
nounced forward scattering, k, = 1.0. 



— ImI]ei,o — Ak/c<; ImSei,i, that describes the elastic 
quasiparticle scattering rate (we have neglected the T^ei,3 
term in ([15]), which vanishes on the Fermi surface, in 
order to highlight the behavior there). In the case of 
a distinct forward scattering potential with = 0.2, it 
is obvious that nodal quasiparticles can be scattered at 
low energies while the scattering of antinodal quasiparti- 
cles is strongly suppressed. In Fig. [21 we show the same 
quantity for a more isotropic scattering potential with 
tv = 1. Here the elastic scattering rate for quasiparticles 
at the Fermi surface shows only a weak anisotropy with 
just a slightly increased scattering rate for nodal quasi- 
particles. With increasing tz^ approaching the isotropic 
limit, the elastic scattering rate at the Fermi surface is 
fully determined by the imaginary part of 5]ei,o while Sei,i 
vanishes. 

In Fig. [3l we compare the energy dependence of the 
nodal and antinodal elastic scattering rates for the two 
different forward scattering parameters discussed above. 
We find for the elastic scattering in the forward scat- 
tering limit {tv = 0.2) a high scattering rate for nodal 
quasiparticles already at low energies compared to the 
gap Aq, while the scattering rate for quasiparticles at the 
antinodal points is strongly suppressed and sets in only 
for higher energies [l^. As previously discussed [18], this 
behavior arises from a cancellation between the "normal" 
So and the "anomalous" Ei channels for co < Ak when 
the scattering is peaked in the forward direction. As ex- 
pected, the scattering becomes more isotropic if the range 
of the impurity potential is decreased {tz = 1) and the 
frequency dependent scattering rates evaluated at differ- 
ent points of the Fermi surface approach each other. In 
this case, the scattering rate for nodal quasiparticles as 
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FIG. 3: (Color online) The elastic scattering rate as a function 
of the quasiparticle energy u for two different values of the 
inverse scattering range n = 0.2 (a) and k = 1 (b). The solid 
Unes show the scattering rate averaged over the Fermi surface, 
while the dashed and dotted lines show the scattering rate for 
quasiparticles at the nodal and antinodal point of the Fermi 
surface respectively. 

well as the Fermi surface averaged scattering rate shows 
a nearly linear increase with frequency, reflecting the d- 
wave density of states. 

B. Elastic ofF-diagonal scattering 

We mention briefly the notion of off-diagonal or pairing 
disorder scattering by out of plane impurities, claimed 
to be responsible in Refs. [H, [H, [H, [U, [H for 
the inhomogeneity in the gap magnitude measured by 
STS, as well as for many of the correlations between var- 
ious STS observables. The scenario described in these 
works is that out of plane dopant impurities, in addition 
to creating a screened Coulomb potential, represent an 
atomic-scale modulation of the BCS pairing interaction 
in the low-energy effective Hamiltonian, via a modula- 
tion of local electronic structure parameters such as t, t' 
or J, or of electron-phonon interaction matrix elements. 
Since the scattering by an induced off-diagonal poten- 
tial ^A(r) is necessarily weak on the scale of the Fermi 
energy, it is again appropriate to use the Born approx- 
imation to describe this scattering. The modulation of 
A(r) observed in experiment [l^, [13, HH is however not 
pointlike, but has a characteristic wavelength of 25-30 A; 
within the picture of Ref. [l2[, this arises because the 
dopant atoms give rise to an atomic scale modulation of 
the BCS pairing potential which causes larger, coherence 



length size fluctuations in Ak(r). This is somewhat dif- 
ficult to model as a quasiparticle scattering potential in 
a disorder averaged calculation (see however [l^), so we 
try in Appendix 1 to estimate the crude effect of order 
parameter modulations on quasiparticle scattering by al- 
lowing the order parameter to be modulated on the 4 
bonds around the impurity, while assuming that it re- 
tains its d-wave character [23. We find, not surprisingly, 
that the scattering rate is largest at the antinode, where 
it rises linearly in uo at zero temperature. 

In an inhomogeneous system like BSCCO-2212, one 
might take ^Ak to be of order Aq to reflect the large 
distribution of observed gap values in the system. The 
overall scale of the scattering parameter F^ defined in 
the Appendix would then vary as Aq, leading to an av- 
eraged local scattering rate of ~ Aqc^;. It is tempting to 
speculate that the measured strong correlation of the STS 
scattering rate coefficient a{v) defined in ([I])-© with the 
local A(r) is due to this effect, but a much more rigorous 
treatment is needed to establish this. 

The results in Appendix 1 give a rough measure of 
the additional scattering rate due to the n or pairing 
disorder component. This would apply, however, only 
in situations where the quasiparticles explore many gap 
"patches" , such that a pairing disorder scenario is valid. 
However, as we argue below, this is the case only at low 
energies, where the STS presents a picture of homoge- 
neous electronic structure. At higher energies, such aver- 
aging is no longer appropriate, and it would be interest- 
ing to consider in more detail the effects of localization 
of quasiparticles within a patch, see Refs. [l2[ and p^ . 
Since a more microscopic approach is not available at 
the present time, and a momentum average will clearly 
give an overall linear-cj contribution to the elastic scat- 
tering rate arising from this term (similar to the ordinary 
Coulomb scattering) , we will ignore its possible relevance 
for the time being and return to it only in the discussion. 



C. Inelastic electron-electron scattering 

Besides the elastic scattering we will also take into ac- 
count the inelastic scattering that arises from the ex- 
change of dynamic spin fluctuations. In the random- 
phase approximation (RPA), the imaginary part of the 
quasiparticle self-energy due to inelastic scattering from 
the on-site Coulomb interaction U can be written as [231 



(16) 
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FIG. 4: (Color online) The sum of the tq and the t\ contri- 
butions to the inelastic scattering rate (which determine the 
scattering near the Fermi surface) as a function of kx and ky 
in the first Brillouin zone for different quasiparticle energies 
cj. In the left column the coupling constant is taken to be 
JJ — 2.20t, which does not produce a well-defined resonant 
mode. In the right column we show results for IJ = 2.36t, 
which produces a resonant mode and an incommensurate spin 
response. From top to bottom cj is chosen to be a; = 0.5Ao, 
uj — Ao and uj = 2.5Ao. The red line shows the Fermi surface. 
The color scale is the same for all six plots and is given below. 
The units of the scattering rate F shown in the color bar are 
Ao. 



In this expression n{(jo) and f{(jo) are the Bose and 
Fermi distribution functions, = y^e^ + and the 
RPA expression for the spin susceptibility x is given by 



X(k,cc;) = 



Xo(k,cj) 
l-/7Xo(k,c^)' 



(17) 



xo(q,^) = 



^k+q^k 



£^k+q£^k 
f Ak+qAk 



/(Ek+q) - /(^k) 



£^k+q£^k 
Ck+q^k + Ak+qAk 
£^k+q£^k 



^ — (£^k+q 

/(i?k+q) - 



+ (-E'k+q ■ 

/(^k+q) ^ 



-£^k) + 
/(^k) 



(-E^k+q + -Ek) 



zO+ J 



where xo is the bare BCS spin susceptibility 



Note that in (p!6|) the coupling constant U is in principle 
different from the effective IJ that enters the denominator 
of the RPA expression Eq. (pT|) for the spin susceptibil- 
ity. Here \J and IJ will be treated as phenomenologi- 
cal parameters. Earlier numerical calculations have used 
U = U = 2.2 to fit microwave and thermal conductiv- 
ity [2^, [32[. Here we will set U = 2.2t and vary U in 
order to see what effects a 7r-resonance produces. 

In Fig. m the sum of the imaginary parts of the tq 
and Ti component of the inelastic scattering rate, which 
determines the inelastic scattering of quasiparticles at 
the Fermi surface, are shown in momentum space for 
the choice U = U = 2.2t (left column) and U = 2.36t 
(right column). Besides the strong scattering peaks at 
the Fermi surface, we find equally strong signals near the 
corner of the Brillouin zone, that are visible for all ener- 
gies 00 and reflect the enhanced spin susceptibility at the 
nesting vectors connecting the two opposite branches of 
the Fermi surface. 

The energy dependence of the scattering rate at the 
nodal and antinodal points on the Fermi surface for U = 
2.2t is now exhibited in Fig. [5](a) and (b). In contrast 
to the elastic rate, the inelastic scattering rate at T = 
is strongly suppressed for low energies and shows an uj^ 
dependence for nodal quasiparticles [H, IH, [sO, IIH, [32[ 
while it vanishes exponentially below Aq for antinodal 
quasiparticles. (Note that in the reference clean nonin- 
teracting quasiparticle system, antinodal quasiparticles 
with cc; < Ao do not exist; the typical or "on-shell" 
lifetime of an antinodal excitation is discussed below.) 
The 00^ dependence of the inelastic scattering for nodal 
quasiparticles can be understood as a combination of the 
Fermi-liquid like uj'^ dependence and the linear uj depen- 
dence of the d-wave density of states. At higher energies, 
on the other hand, the broadening of the antinodal states 
in the case of weaker interactions (Fig. [5] (a) and (b)) be- 
comes somewhat larger than that of the nodal states, due 
to the crossover to the normal state band structure. In 
Fig. [5] (b) we have compared the energy dependence at 
the node and antinode with the Fermi surface average, 
which will be useful to us below. It is worthwhile not- 
ing at this stage that the full Brillouin zone average of 
rinei(k, cj) (not shown) does not differ qualitatively from 
the Fermi surface average. 

Within the generalized RPA spin fluctuation approach 
to the Hubbard model, it has been known for some time 
that if U is adjusted to tune the system closer to the 
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FIG. 5: (Color online) The inelastic scattering rate rinei(k, u) 
for k on the Fermi surface as a function of the quasiparticle 
energy a; for different values of interaction parameter U — 
2.20t (a,b); 2.36t (c,d); and 2.56t (e,f). Left panels: scattering 
rate at different k ranging in equal increments along Fermi 
surface from node (blue) to antinode (red). Right panels: The 
solid lines show the scattering rate averaged over the Fermi 
surface, while the dashed and dotted lines show the scattering 
rate for quasiparticles at the nodal and antinodal point of the 
Fermi surface respectively. A constant broadening of 0.1 Ao 
has been added to regularize the numerical calculation. 



antiferromagnetic phase transition, a resonant S = 1 col- 
lective mode is pulled down below the particle-hole con- 
tinuum in the dynamical spin susceptibility x" {^.i ^) near 
q = (tt, tt). This mode has been identified as the sharp 
spectral feature observed in neutron scattering measure- 
ments at this wavevector [33| . The RPA represents one of 
several approaches which have been proposed to describe 
the TT-resonance which is observed in the cuprates [U- 
Specific applications to the neutron response of different 
cuprate materials have been given in Ref. [35[. 

Here we first discuss the basic kinematics of the collec- 
tive mode and consider the effect on the fermionic one- 
particle self-energy within our model. This has been con- 
sidered already using approximate models for the suscep- 
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FIG. 6: (Color online) Left: imaginary part of dynamical 
susceptibility x"{^.^^) vs. uo for the three values of U consid- 
ered in the text. Right: false color plot of x"{^.^^) for a cut 
through the Brillouin zone along (110). 



tibility [36|, |37|, |38|, l39| for purposes of comparing features 
of the neutron and ARPES spectra. Within the RPA 
form Eq. [T71 a resonance in the real part of the suscepti- 
bility may occur when 



U = l/X0(Q,^res). 



(18) 



For our band parameters, and if we require Q = (7r,7r), 
the resonant frequency occurs at uOres = 0.9Ao for U = 
2.56. More generally, for Q sufficiently close to (7r,7r), it 
may be shown that XoiQ^^res) vanishes at T = 0. In 
this case there is a sharp collective mode contribution of 
the form 



x" 

Ares 



(Q,w) 



(19) 



to the imaginary part of the spin susceptibility near Q. 
In Figure [6l we show how, with increasing interaction 
strength, the energy position of the collective feature in 
X^^ shifts downward and the intensity increases. These as- 
pects of the spin fluctuation spectrum are then directly 
reflected in Fig. [5l Many of the features of the experi- 
mental neutron spectrum on optimally doped YBC0[4l| 
and BSCCOjiil, including the hourglass-like shape of the 
neutron response in q, a; space (Fig. [6] (b)), as well a 
characteristic rotation of the pattern of incommensurate 
q-peaks as the energy is tuned through the resonance are 
known to be captured by this type of theory [s^, EO] • 

To see how the collective mode affects the quasiparti- 
cle scattering rate, we focus on the resonant part of x^' . 
The contribution to the effective interaction l^res(q, ^) is 
^U'^Xresi^^^)^ giving rise to an additional contribution 
to the total scattering rate 



Tres{K(^) = -l 5^ ImKes(k-k^cJ + s^kO(20) 



(/(5^kO + ^(^ + 5^kO) 



1 - 5- 



At the node k = kAr, the expressions simplify consid- 
erably and we find at T = that 



8 



4U 



^S{l-Uxo{^-k',uj-Ei,,)) 



(21) 



T he integr and contributes provided co = cOres + 
ye^T+Aj^, with = k — Q. Since for our band pa- 
rameters k' = kAT — (tt, tt) is far from the Fermi surface, 
the contribution to T from the magnetic resonance mode 
win arise at frequencies uo = ujres + kk'l larger than cOres 
itself by an amount which is of order a fraction of the 
Fermi energy, much larger than Aq. 

For the antinodal points, or rather the nearby "hot 
spots" ]<.HS where the Fermi surface crosses the an- 
tiferromagnetic zone boundary, the wave vector k' = 
^HS — (tt, tt) is also on the Fermi surface. However the 
off-diagonal parts of the coherence factors must now also 
be included, so the resonant contribution becomes 



-L res 



417 



x0{lo - Ek,) 



^(5(l-C7xo(k-k',a;-Ek')) 



(22) 



where /S^hs ^ Aq is the value of the order parameter 
at the hot spot k//^. Thus there is a contribution to 
the resonant part of the scattering rate at the hot spot 
when uo uOres + Aq. This additional scattering has 
been claimed to be responsible for part of the "peak-dip- 
hump" structure seen in ARPES near the antinode at 
low temperatures in the superconducting state [H, . 
However, since uOres is also of order Aq, it does not appear 
as though the resonant electron-electron scattering can 
by itself play a significant role in the electronic scattering 
rate for energies at or below the coherence peak energy 
Ao. 

Nevertheless, the enhanced inelastic scattering will 
play a role in measured STS properties. The tuning of 
the system closer to the antiferromagnetic instability has 
the effect of enhancing the overall scattering rate, such 
that the states near the nodes now decay more rapidly. 
In Figure [5jc)-(f), we show the evolution of the energy- 
dependent scattering rate for states along the Fermi sur- 
face for slightly larger U . The contribution of the col- 
lective mode is easiest to discern in (e)-(f), for the case 
IJ = 2.56. The scattering in the midgap range near the 
node has risen considerably, although near the antinode 
it is still exponentially suppressed at low energies. The 
collective mode results in a contribution which rises near 
uo :^ cJres + Ao ~ I.QAq at the antinode, as expected from 
the above discussion, which then peaks at about 2.5Ao. 
There is no obvious resonant contribution at higher en- 
ergies in the nodal direction, but its existence can be 
deduced by considering the evolution of the antinodal 
peak as one moves away from the antinode: as seen in 
the figure, it moves to higher energies as expected and 
broadens due to the additional decay channels. Above 



an energy of order 1.5Ao, the antinodal scattering rate 
is seen to become much larger than the nodal scattering 
rate. This enhancement is directly related to the open- 
ing of the superconducting gap and the creation of the 
resonant mode. 

The intermediate case U = 2.36 (Fig. [5Kc)-(d)) also 
shows clear evidence of enhanced interactions, and ves- 
tiges of the mode features. The effect on Finei at low 
energies is primarily quantitative: the inelastic scatter- 
ing near the node becomes significant in the midgap 
range compared to the weaker case (a)-(b), and is al- 
ready a fraction of Aq near the node at the gap frequency 
cc; = Aq. This will play a role in our discussion of STS ex- 
periments, because the extracted quasiparticle scattering 
rate is of this order for optimally doped samples. There 
is further independent evidence in the work of Fauque 
et al. [42] to suggest that the intermediate interaction 
strength case we considered here is most relevant for ex- 
periments on optimally doped BSCCO. This is because 
in Ref. [42] , the collective mode energy is measured to be 
42 meV, while the average gap is extracted as 35 meV. 
This ratio of about 1.2 is achieved in our model for inter- 
action strength IJ = 2.36; changing this parameter moves 
the resonance up or down relative to the gap edge, as seen 
in Fig. [nja). 



D. Total scattering rate 

To include both elastic and inelastic scattering ef- 
fects, we neglect interference processes between electron- 
electron collisions and impurity scattering entirely, and 
approximate the total scattering rate by 



inel- 



(23) 



In Fig. [7J the total scattering rate for quasiparticles at 
the nodal point of the Fermi surface (dashed) is compared 
to the scattering rate for quasiparticles at the antinodal 
point of the Fermi surface (dotted) as well as the average 
scattering rate including only the quasiparticles at the 
Fermi surface. 



■ avg 



(cj) = (rtot(k,ct;)) 



(24) 



where i^ps is a Fermi surface average. The total scat- 
tering rate is seen in Fig. [3 (a)-(b) to be dominated for 
uj < 2Ao by the linear energy dependence of the elas- 
tic part, while at higher energies uo > 3Ao it refiects the 
quasilinear energy dependence of the inelastic scattering 
rate in the classical regime. The slope for low energies 
and the slope for high energies are not related to each 
other and depend on the individual parameters used in 
the two models. 

For completeness, in Fig. [3(b) the total scattering rate 
averaged over the Fermi surface is shown together with 
its two averaged contributions from elastic and inelastic 
scattering processes. It shows qualitatively the same en- 
ergy dependence as the Brillouin zone averaged rate (not 
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FIG. 7: (Color online) The elastic, inelastic, and total scatter- 
ing rates as a function of the quasiparticle energy uj for several 
assumptions about the strength of each scattering channel. In 
the left panels, the total rate Vtot at the antinode (node 
Icat) is plotted as a dotted (dashed) line. The average over 
the Fermi surface of rtot(k,a;) is plotted as the solid line. In 
the right panels, the contribution of the elastic (blue) and 
inelastic (orange) to the total (magenta) Fermi surface av- 
eraged rates is exhibited separately, (a)-(b) are for the case 
IJ = 2.36t with Vei,N(Tc) = O.IAq. (c)-(d) Same as (a)-(b) 
but with rei,N(Tc) = 0.5Ao. (e)-(f) Same as (a)-(b) but with 



shown) but has a somewhat larger overall magnitude. 
This is due to the fact that the Brillouin zone average 
includes regions of the zone away from the Fermi surface 
with small scattering rates. 



IV. LOCAL GREEN'S FUNCTION 

In this section we calculate the exact local Green's 
function including impurity and spin-fluctuation scatter- 
ing, and exhibit the LDOS for varying amounts of dis- 
order and different interaction strengths. Our first goal 



is to see if the use of a local self-energy, which we define 
in an ad hoc way, can provide a good description of the 
exact LDOS calculated with the full self-energies. Sec- 
ondly, we compare these approximations to the linear-cj 
ansatz, given by Eqs.[T][2l employed by the Cornell group 
to fit their data [8]. We would like to see if an approxi- 
mate local self-energy may be constructed from the local 
Green's function, and to what extent the true self-energy 
yields an LDOS similar to that obtained from the linear- 
00 scattering rate. 

The local Green's function G(r, r, cj) in the presence of 
elastic and inelastic scattering is given, within a region of 
size of order ^, as the Fourier transform of the momentum 
dependent Green's function G(k, uo) 



(25) 



Here the full momentum dependent Green's function is 
calculated from 



o^tq + 6kr3 + AkTi 



(26) 



and a), ek and Ak are the renormalized quasiparticle en- 
ergy, band structure and gap in the presence of the elas- 
tic and inelastic self-energy calculated in the previous 
section. The LDOS is then simply given by 

iV(r, = -^Tr [{to + ts)G{v, r; uj)] . (27) 

The results of these calculations including the full self- 
energies are shown as the solid lines in Fig. [8] for (a) pure 
elastic scattering, (b) pure inelastic scattering and (c) for 
the total scattering rate. The green curve shows for com- 
parison the the LDOS for a clean d-wave superconductor. 
The high slope of the elastic scattering rate leads to a 
linear energy dependence of the total rate up to a value 
of Ao and a strong suppression of the coherence peaks. 
The inelastic scattering rate, which is strongly suppressed 
for low energies, begins to significantly affect the LDOS 
only at energies a; > Aq for the slightly off-resonant case 
U = 2.36t case shown, and has qualitatively no influence 
on the low energy spectrum. In the same figure panels, 
we compare these results to calculations with the ad hoc 
"local" scattering rate Tavgi^) which we simply add as 
an imaginary part to the quasiparticle energy, 

G [r, r, + iTavg (cj)] = ^ G [k, + iTavg (oo)]. (28) 

k 

As can been seen, this approximation (dashed lines) leads 
to a slightly lower spectral weight for low energies but 
exhibits also a quasilinear energy dependence for the local 
density of states up to an energy of Aq. Finally we fit 
the effective momentum averaged scattering rate with the 
linear energy dependence of Eq. [2] and we find for low 
energies excellent agreement with the approximation of 
a momentum averaged scattering rate (orange). 
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FIG. 8: (Color online) The local density of states (LDOS) 
A^totai(^) normalized to the Fermi level density of states A^o 
for a d-wave superconductor with scattering processes of var- 
ious types included. In all panels, green line shows LDOS for 
pure system, blue (a) gives the LDOS for pure elastic scat- 
tering with Fei (Tc) = O.lAo, red curve (b) gives pure inelas- 
tic scattering with U = 2.36, and magenta (c) includes both 
elastic and inelastic processes. Solid lines are calculated with 
the full momentum dependent self-energy, while the dashed 
curves used an effective momentum averaged scattering rate 
defined in (|24p . In (d) the LDOS in the presence of impurity 
and spin fluctuation scattering, with parameters from (a) and 
(b), is compared to the best fit Niin(oo) (orange) using a lin- 
ear scattering rate of the form given in Eqs. {[T]|2]) as extracted 
from STS data in Ref. 8. Figs, (e) and (f) show the LDOS 
for a dirtier system with Fei (Tc) = 0.5Ao (line types same as 
in (a)-(d)). 



Thus we conclude that for energies co < Aq^ there is 
a quasilinear increase of T arising from elastic scattering 
processes and in calculating the LDOS it can be mod- 
eled by the linear form Eq. ([2]). On the other hand, the 
STS local scattering rate is also enhanced by inelastic 
scattering processes when a; > Aq. Despite the fact that 
this destroys even the approximate linearity of the total 
scattering rate (see Fig. [7|), the momentum average in- 



herent in the local measurement implies that it is very 
difficult to distinguish the true functional form of F in 
the impurity plus spin fluctuation model from the phe- 
nomenological linear ansatz of Ref. \§\ . We return to this 
point in Section \V\ below. 

The results in Fig. [8] are qualitatively similar to earlier 
attempts to fit STS data by assuming a model self-energy 
dressing a BCS-like d-wave Green's function. In particu- 
lar, Hoogenboom et al. [i^ and de Castro et al. Q com- 
pared data on BSCCO to self-energy models with con- 
stant scattering rates, a marginal Fermi liquid model [45[, 
and the model of Ref. [s^ describing a coupling to a 
phenomenological collective spin mode. They concluded 
that details of the peak-dip hump structure observed in 
some spectra required the collective mode, and identified 
features at sum and difference energies of the gap and 
collective mode frequency. Here we have attempted no 
detailed fits of STS data, as we are primarily interested 
in exploring the general notion of a local scattering rate. 
We note that Fig.[8]shows, as in Refs. 0,113, the broad- 
ening of the van Hove singularity by inelastic scattering, 
and that out-of-plane elastic scattering may play a role 
as well. Peak-dip-hump structures similar to experiment 
are seen in some cases in Fig. [51 but more investigation 
is needed to determine whether these aspects of the ex- 
perimental spectra correspond to scattering effects of the 
type considered here. 

V. COMPARISON WITH ARPES SPECTRAL 
FUNCTION 

The model of the quasiparticle scattering rate put for- 
ward above has implications for other quantities besides 
the STS conductance, and we may learn something from 
the comparison. In particular, as mentioned in the in- 
troduction, ARPES has reported superconducting state 
relaxation rates which are many times larger than those 
measured by STS. A naive application of these ARPES 
rates in the superconducting state would lead to even 
broader LDOS spectra than those in Fig. [8](e)-(f), in 
contradiction to experiment. The ARPES signal is given 
in the sudden approximation by 

/(k,t^) = |MpA(k,c^)/(c^), (29) 

where M is an ARPES matrix element, which we take to 
be constant, 74(k,a;) = — (l/27r)ImTr (tq + T3)G(k, o;), 
and / is the Fermi function. Here we calculate the 
full electron spectral function ACk^co) within our sim- 
ple model, for different points on the Fermi surface and 
for different scattering rates, with the aim of describing 
some qualitative phenomena and comparing with STS. 
We note that the Fermi surface itself is, in the presence 
of interactions and scattering, not identical with the non- 
interacting pure one, but rather is the solution in k space 
of the equation 

Re ik = ek + Re {^f + E^^^^ ) = 0. (30) 
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FIG. 9: (Color online) Left: spectral function A(k,uj) for k on 
the Fermi surface at T = with elastic scattering parameters 
= 1 and (a) rei,N{Tc) = O.lAo and (c) Fei^NiTc) =_Ao. 
and with the inelastic scattering rate calculated with U = 
2.36t. Color ranges in equal increments from blue (node) to 
red (antinode). Only four k points are shown in (a) due to 
narrowness of peaks. Right: spectral gap as determined by 
one-half the energy difference between the two maxima of 
A(k,cj) for (b) rei,N{Tc) = 0.1 Ao and (d) = Ao, respectively. 
Red curve in both figures indicates pure d-wave gap in the 
absence of scattering. 



Here the chemical potential is determined to maintain 
a given filling. We have verified that this does not af- 
fect the broadening and in the calculation of the spec- 
tral functions shown in Fig. [9] the renormalization of the 
Fermi surface by Re S was neglected. 

In Fig. [9] (a),(c), we show the evolution of the low- 
temperature spectral function along the Fermi surface, 
beginning with the (blue) peak centered at = as 
expected, and ending with two clear (red) peaks at the 
antinode. Close to the node, some spectra still have their 
maxima at the Fermi surface due to scattering effects. 
In the pseudogap state of the cuprates, a similar phe- 
nomenon is observed: the maximum of the spectral func- 
tion remains at the Fermi surface over some range of k 
values centered on the nodal point; the pseudogap itself 
is visible as a double-peak structure only some distance 
away from the nodal point. The range of k where the sin- 
gle maximum is at the Fermi energy is called the Fermi 
"arc" , and it evolves continuously into the full Fermi sur- 
face in the normal state as the pseudogap disappears. It 
has been pointed out by several authors [El, @, 0] that the 
arc phenomenon can be explained trivially by appealing 



to an energy and/or temperature-dependent scattering 
rate, which for sufficiently large scattering broadens the 
two-peak spectral function one would expect in the pres- 
ence of a spectral gap into a single peak centered at the 
Fermi level. A similar phenomenon occurs here in the su- 
perconducting state, as shown in Fig.[9](b),(d), where the 
effective gap, determined from the position of the spec- 
tral function peak, is plotted along the Fermi surface, 
exhibiting a finite range of k points where it vanishes. 
For a constant scattering rate F, the criterion determin- 
ing the angular position of the end of the arc is F = aAk, 
where a is a constant of order unity, equal to ^/3 for a 
d-wave superconductor with circular Fermi surface and 
cos2(/) order parameter 0. Note the effective gap deter- 
mined by ARPES in this manner does not correspond 
to the renormalized order parameter in the theory, as 
pointed out by Sensarma et al. [i^. 

Since the scattering rate is temperature and frequency 
dependent, the criterion for the position of the spec- 
tral peaks changes as T is increased, resulting in a 
temperature-dependent change in the "arc length", as 
shown in Fig. [TOl This is a very weak temperature de- 
pendence until quite close to Tc, since it is driven by the 
T-dependence of the gap in the theory [19], and above 
Tc the arc length in the current model is fixed because 
the gap goes to zero. A theory including a pseudogap 
in the normal metallic dispersion would give rise to a 
T — dependent (linear-T for a marginal Fermi liquid-like 
scattering rate[45|) arc length [H, H, 0]. Here we have 
shown that, at low T in the superconducting state, Fermi 
arcs are also possible, but only in a dirty system. Thus 
the recent observation by Kanigel et al. [47| of the col- 
lapse of the arc length at low T in the superconduct- 
ing state places constraints on the elastic scattering rate. 
Since no arc is observed at low T in optimally doped 
samples, the scattering rate must be a tiny fraction of 
Ao, inconsistent with ARPES determinations of such 
rates at low T, but consistent with the STS result Q 
of T{u; Ao) ~ 2 meV near optimal doping. On the 
other hand, STS has shown that the linear fits to the 
scattering rate yield a tenfold increase in this rate for 
highly underdoped samples. In this case, F becomes an 
appreciable fraction of Ao for subgap energies. The im- 
plication is therefore clear that ARPES should eventually 
observe arcs in the superconducting state as the system 
is underdoped. 

Finally, it is worth observing that in Fig. [9] (a), (c) the 
EDC width at the antinodal points is significantly larger 
than at the antinode in the clean case shown in[9](a)), 
despite the fact that the antinodal scattering rate does 
not exceed the nodal one until energies several times Ao. 
This is primarily because the EDC's are dominated by 
quasiparticles "on the mass shell" co = Ak- Thus, while 
the nodal scattering rate is actually larger than the antin- 
odal scattering rate for u = Ao, as shown in Fig. [71 
the scattering rates which broaden ACk^co) are associ- 
ated with on-shell co = A^ energies. In this case, the 
antinodal scattering rate refers to quasiparticles of en- 
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FIG. 10: (Color online) Length of "Fermi arc" defined in text 
at T = in superconducting state vs. temperature T/Tc. 
Elastic scattering rate parameter Tei^N — Ao 



ergy ~ Aq, while the nodal scattering rate to those of 
energy a; ~ 0. In addition, if interactions are sufficiently 
strong, the inelastic scattering at the antinode will even- 
tually contribute, leading to a pronounced asymmetry 
of the spectral function in the form of a large tail on the 
high binding energy side in ARPES. This may contribute 
significantly to measured EDC widths if electron-electron 
scattering is sufficiently strong. 



VI. EFFECT OF SAMPLE INHOMOGENEITY 
A. Low energy homogeneity 

Thus far we have considered a homogeneous d-wave 
superconductor with inelastic scattering and microscopic 
disorder, which we treated by averaging to obtain a 
translationally invariant self-energy characterized by mo- 
mentum k. We saw in the previous section that at 
a given quasiparticle energy, calculating a "local self- 
energy" S(k, cj) provided a reasonable description of 
the LDOS broadening. This assumes, as discussed in the 
introduction, that the system is self-averaging within a 
region with a single gap magnitude Ak- However, the 
BSCCO-2212 system is known to be inhomogeneous at 
the nanoscale, in a somewhat mysterious way. Spectra 
measured near the typical gap edge of a given sample are 
extremely inhomogeneous, with coherence peak positions 
varying up to a factor of 2-3 within a typical experimental 
field of view. On the other hand, below some typical en- 
ergy of order 30 meV, near- homogeneity is recovered [IJ] • 

Let us consider the STM-determined low-energy lo- 
cal scattering rates. These are quite small relative to 
rates determined by ARPES, but increase with increas- 
ing energy. It is our hypothesis that quasiparticles at 
low energy scatter so seldom that they explore many gap 
"patches", and thus do not carry knowledge of the lo- 



FIG. 11: (Color online) Quasiparticle group velocity for band 
structure given in ()12p for a) normal state; b) superconducting 
state. Length of arrows or grayscale indicate size of local 
quasiparticle speed. 



cal gap value when collected by the STM tip. As the 
energy increases, eventually their mean free path drops 
until they are localized within a single patch. Conduc- 
tance spectra taken at these higher energies will therefore 
reflect the local gap at the position of the STM tip. It is 
significant (a) that the energy where inhomogeneity ap- 
pears is of order Aq in Bi-2212, and b) that the doping 
dependence of this energy scale is rather weak. 

To calculate the energy-dependent mean free path of 
our system, we will need the velocity of a typical Bogoli- 
ubov quasiparticle v^p at different points near the Fermi 
surface at T = 0. The group velocity for a quasipar- 
ticle of momentum k in a clean BCS superconductor is 
= V£^k- Here however we are accounting for scatter- 
ing of these quasiparticles by impurities and collisions; 
therefore there is an effective spread in momentum of a 
typical quasiparticle due to its lifetime broadening. We 
therefore define the speed of a typical quasiparticle near 
the Fermi surface with momentum k = {k±^k ||) and 
energy at the gap edge to be 



(31) 



At the node, the spectral function is (neglecting real 
parts of self-energies) ACkN^co) :^ (r/7r)/(^^ + F^), so 
using (|3T]) we find ^^^(kAr) = vfO^n)- At the antinode, 
the spectral function is approximately 

F 2Ao (Ao + 



^4 + 4A2F2 



leading to 



.gp(kA) ^ \ H- vfO^a), 

^0 



(32) 



(33) 



where we have assumed F <C Aq. Thus the "typical" 
antinodal velocity is generally much smaller than the nor- 
mal state Fermi velocity. Using the STM determined 
value of F = 2 meV at optimal doping, Vqp should be 
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FIG. 12: (Color online) Typical quasiparticle group velocity 
for band structure given in ([12]) vs. position on Fermi sur- 
face using expression (|3T]) . The inelastic scattering is calcu- 
lated for U = 2.36t while the elastic scattering corresponds to 
riv(Tc) = 0.05Ao (a) and F = 0.2Ao (b). 



suppressed by a factor of 3-4 relative to the Fermi veloc- 
ity. This is confirmed by a full evaluation of (|3T]) using 
the self energy given by Eq. [23l 

We now define the on-shell quasiparticle mean free 
path 



4 = 



2F(c^ = Ak)' 



(34) 



where F may represent the elastic, inelastic, or total one- 
particle scattering rate. For extremely dirty systems, mo- 
mentum is not a good quantum number, and the on-shell 
quantities may not provide a good representation of the 
mean free path. However for the doped cuprate sam- 
ples of interest it appears to us that this is a reasonable 
estimate. The on-shell scattering rate is precisely that 
defined for the cj— dependent momenta k corresponding 
to the points of the contours of constant quasiparticle 
energy which lie on the Fermi surface, i.e. the so- 
called "banana tips" defined within the octet model [i^. 
In Fig. [131 we plot both the scattering rates and mean 
free paths defined in this way. These are shown both 
as a function of uo and, equivalently, of k, for the set of 
parameters corresponding to a clean system with inelas- 
tic scattering rate parameter U = 2.36t consistent with 
neutron data, as in Fig. [7] (a) and (b) and Fig. [8] (a) and 
(b). In Fig. [13] (a) and (b), we see that the effective on- 
shell scattering rate is indeed quasi-linear in energy over 
the range < u < Aq^ resulting in a momentum depen- 
dence similar to the d-wave gap itself. The parameters 
chosen to approximately reproduce the neutron scatter- 
ing behavior in YBCO ( Fig. [3 (a)-(b)) appear to give a 
scattering rate scale of approximately O.IAq, of order the 
correct scale of 2meV determined by STS for BSCCO. 

The mean free path in Fig. [13] is seen to fall rapidly 
with increasing energy, and becomes of order the patch 
size ~ 10a at an energy around the gap edge. We will des- 
ignate this energy cuioc since it is only above this energy 
where true local behavior characteristic of a particular 
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FIG. 13: (Color online) (a) on-shell scattering rate rtot(k, co = 
Ak) as a function of k from the node to the antinode along 
the Fermi surface parameterized from -1 to 1, for scattering 
parameters elastic scattering parameters = 1, Fei (Tc) = 
0.1 Ao and U = 2.36t. Blue: elastic scattering; red: inelastic; 
magenta: total scattering rate, (b) same quantity plotted vs. 
(jj. (c,d) Mean free path i = t'gp/(2rtot(k, u; = Ak)) vs. k and 
CO as in (a,b). 



gap patch order parameter is measured. This is reminis- 
cent of the behavior identified by STS on BSCCO-2212 
samples: below a critical energy uoioc of order or slightly 
less than Aq, spectra are homogeneous. The famous in- 
homogeneity measured by STS in this system appears 
only above this energy, and in particular at and around 
the average gap energy (and above for underdoped sam- 
ples). The energy coioc above which the sample inhomo- 
geneity is sensed by the STS measurement is therefore 
where spectra measured in different locations begin to 
differ from one another, and must therefore necessarily 
be associated with a change in slope N{u;) at coioc Such 
a kink in typical spectra where homo gene ity is lost has 
indeed been observed in STS data [H,]!^. 

We note also the relative insensitivity of the observed 
kink energy to doping pjsO], despite rapid changes in both 
the scattering rate and gap size over the measured range 
of p ^ 0.08 — 0.2. Within our picture, this can be un- 
derstood at least in part by the increase of the antinodal 
typical quasiparticle velocity with the scattering rate as 
shown in Eq. ([33]) . Changes in Fermi surface shape may 
also play a role in keeping the mean free path at this 
energy roughly doping independent. 

Finally, we comment on the significance of these results 
for the quasiparticle interference (QPI) patterns observed 
in Fourier transform STS on these systems [4^. It is im- 
portant to recall that the QPI patterms arise due to the 
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random static potential from impurities. Static disorder 
can enhance noise which may then swamp certain Fourier 
transform q-peaks[5l|, but it cannot suppress or broaden 
those peaks in q-space. Inelastic scattering however will 
dephase quasiparticle wavefunctions, and it is to be ex- 
pected that a given FT-STS q-peak will be lost when the 
inelastic scattering rate becomes large enough such that 
the de Broglie wavelength of a quasiparticle at the Fermi 
level is smeared by a significant fraction. The theoretical 
question of how such quasiparticle interference patterns 
are destroyed is an interesting one, which we are currently 
studying in more detail. At present, however, our results 
suggest that a qualitative increase in this smearing oc- 
curs in the neighborhood of cj* ~ Aq due to the rapid 
rise of the inelastic scattering. It is not clear from our 
analysis that the two energy scales coioc where the system 
inhomogeneity becomes manifest and the scale a;* where 
QPI patterns are destroyed are the same scale, but they 
are of the same order, and experimentally appear to be 
quite close (5Q|. 

B. STS and ARPES quasiparticle relaxation rates 

The size of the ARPES laser spot on the sample sur- 
face is many times larger than a typical "gap patch" size 
in BSCCO-2212. This fact has been frequently pointed 
out, but its consequences for extracted quasiparticle life- 
times has not been explored to our knowledge. We noted 
in the introduction that ARPES-extracted lifetimes for 
both nodal and antinodal quasiparticles are much larger 
than those determined recently by STS [8]. This fact 
can be plausibly explained by taking our current knowl- 
edge of STS lifetimes and averaging the spectral func- 
tion over the distribution of gaps P(A) found near the 
BSCCO-2212 surface. A similar analysis was performed 
by Fauque et al.[42] to explain the large width of the neu- 
tron resonance peak in BSCCO-2212, assuming that the 
same distribution characterized the bulk of the sample 
probed by neutrons. The gap distribution in STS has 
been found to have a form roughly equal to a Gaussian 
centered at the average gap A and having a width a ap- 
proximately equal to 0.2 A, more or less independent of 
doping [8]. In Alldredge et al., the STS scattering rate at 
the gap energy was extracted as ~ acj, as described in the 
introduction. The average a for an optimally doped sam- 
ple was {a) c:^ 0.04, giving a scattering rate at cj = A 
40meV of - 2meV, or 0.05A. 

The spectral function at the antinode where this gap 
value is achieved is shown in Fig.[T4las the dashed line, a 
narrow Lorentzian centered at Aq. But ARPES averages 
over many such lines, each centered at a different Aq and 
broadened by a different scattering rate, which depends 
on Ao in the above analysis. The solid line in Fig. [14] is 
then given by 

AARPEs(k,c^) - [ P(Ao)A(k,cj; Ao), (35) 
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FIG. 14: (Color online) T=0 spectral function at antinode 
for homogeneous d-wave superconductor with scattering rate 
of Pel (Tc) = 0.05Ao (dashed lines). The solid line is the 
same spectral function convolved with a Gaussian distribution 
of gap values P(Ao) taken from experiment [sl]. Full width 
at half-maximum of the solid curve is about 5x that of the 
dashed curve. Note no instrumental energy broadening has 
been added to this curve. 



where 74(k, cj; Aq) is the spectral function calculated in 
Section [Vl Note to make this plot inelastic scattering was 
neglected in order to make the computation time prac- 
ticable. However our investigations showed that the dis- 
tribution of center positions (gap values) is in any case 
much more important for the ultimate measured EDO 
width than the distribution of (A-dependent) scatter- 
ing rates. This effective spectral function is roughly five 
times wider than the "intrinsic" spectral function which 
would be measured locally in a small region where the 
local gap A(r) is A. Together with the instrumental res- 
olution of ARPES near the antinode of roughly lOmeV, 
which should be convolved with Eq. (|35|). this gives an 
effective antinodal EDO full width at half maximum of ~ 
25 meV. This is very similar to the full widths of antin- 
odal ARPES EDC's currently reported[^. It therefore 
appears to us that the inhomogeneity effectively prevents 
ARPES (and planar tunnelling measurements, which suf- 
fer from the same problem) from measuring intrinsic life- 
times in the superconducting state of the BSCCO family 
of materials. 



VII. CONCLUSIONS 

In this paper we have examined a series of questions 
arising from recent STS measurements of lifetimes in 
BSCCO extracted from fits of the conductance to a BCS 
d-wave form with broadening, Eqs. ([T]) and ([2]). We have 
argued that an STS measurement probes local d-wave 
quasiparticle states averaged over a mean free path, with 
the exception of unitary bound states where quasiparti- 
cles are trapped on an atomic scale. We have discussed 
different mechanisms for the scattering of quasiparticles 
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in the superconducting state, and presented the results 
of model calculations within a BCS framework. Among 
these scattering processes are elastic quasiforward weak 
scattering by out of plane dopants, elastic pairing disor- 
der scattering by the same dopants, and elastic unitary 
isotropic scattering by in-plane defects. We also consid- 
ered the effects of inelastic scattering by spin fluctuations 
treated within a generalized random phase approxima- 
tion, for different values of the quasiparticle interaction 
parameter, corresponding to a dynamical spin suscepti- 
bility at TT, TT with a resonant mode contribution in the 
(i-wave state. 

Within this model, the effective scattering rate rises 
with energy, and is generally linear at low energy due to 
weak elastic processes. If interactions are weak, quasi- 
particles over the whole "low-energy" range uo < Aq are 
scattered primarily by impurities. If the system is suf- 
ficiently clean, or if quasiparticle interactions are suffi- 
ciently strong such as to induce a resonant mode, one may 
have a situation where the elastic processes dominate 
only up to a subgap crossover energy ^ Aq. In the case 
with resonant mode considered in Fig.[5je)-(f), for exam- 
ple, this energy was of order O.SAq and the inelastic scat- 
tering rate for near-nodal states rose rapidly above this 
energy; weaker interactions move this crossover above the 
gap edge. The energy dependence of the total scattering 
rate for any of these situations is never strictly linear as a 
function of uo. However, we found that the model LDOS 
is well approximated by the imaginary part of a d-wave 
Green's function with a linear self-energy. This was un- 
derstood in terms of the effective on-shell total scattering 
rate, which indeed appears to be quasilinear for reason- 
able parameters. Therefore, while the STS extraction of 
local lifetimes is not particularly sensitive to the details 
of scattering of the low energy states, it is roughly cor- 
rect and provides important insight regarding the doping 
dependence of local scattering. 

In this work we have focussed primarily on optimal 
doping, where a BCS weak-coupling approach may be 
expected to work well. We have therefore not addressed 
this striking doping dependence of the STS scattering 
rate, where increases of a factor of ten or more relative 
to optimal doping were reported in strongly underdoped 
samples [H. Within the present theory, such effects might 
be captured if the effective U /t were to increase with un- 
derdoping, but this is difficult to describe in a systematic 
manner in such a framework, as is well known. In ad- 
dition, we have not addressed the measured correlation 
of the scattering coefficient a with local gap size Aq. In 
fact, within the description of elastic potential scattering 
discussed above, an anticorrelation would be obtained, 
in contradiction to experiment. On the other hand, pair 
disorder scattering, which we have also discussed briefly, 
appears to provide a qualitatively correct description of 
this empirical fact, as well as a quasilinear in uj local 
scattering rate. While we have not provided a consistent 
microscopic description of this physics, it appears to us 
to be worth exploring further. 



As energy increases, the rising quasiparticle scattering 
rate leads to a falling mean free path. Our estimates 
show that at low energy, quasiparticles explore (and self- 
average over) a large area, providing an explanation of 
why STS spectra are homogeneous at low energies. At 
some critical energy uoioc of order Aq, the mean free path 
becomes of order the gap patch size, and STS then (for 
^ > ^loc) measures spectra broadened by a local self- 
energy characteristic of a single gap. The fact that a 
change in slope in conductance spectra is observed at 
nearly the same energy in all patches appears to suggest 
that the scattering rate rises fairly rapidly in this energy 
range. Above this energy, conductance spectra are no 
longer homogeneous since self- averaging over many gap 
patches no longer takes place. It is furthermore to be 
expected-although we have not shown this directly-that 
dispersive quasiparticle interference patterns observed in 
Fourier transform STS will disappear at roughly this en- 
ergy as well, since their existence depends on the as- 
sumption of a quasiparticle with definite energy scatter- 
ing many times from a disorder potential. 

To compare with ARPES, we have used the same 
model to calculate the spectral function A{\^^uo)^ and 
shown that a "Fermi arc" — where the spectral peak re- 
mains at the Fermi level away from the node in mo- 
mentum space — may exist in the superconductin state. 
We argued that while current experiments on optimally 
doped samples do not see this feature, it should be- 
come visible in underdoped samples as scattering rates in- 
crease. Actual ARPES EDC linewidths are much broader 
in energy than those found within our model of a homoge- 
neous d-wave superconductor. We have therefore argued 
here that this is due in large part to shifts of the spec- 
tral function caused by gap inhomogeneity, at least in the 
BSCCO family of cuprates, and estimated the effective 
EDC width which should be measured in such a system. 
Additional effects which add to this large enhancement 
of EDC linewidths include ARPES instrumental broad- 
ening and possible bilayer effects near the antinode. The 
current theory appears to account for the much larger 
ARPES EDC widths in the superconducting state com- 
pared to STS. 

In this paper we have concentrated primarily on quali- 
tative physics, leaving open questions of quantitative fits 
to the ARPES and STS spectra. We hope to address 
these in a future publication. 
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VIII. APPENDIX 1: MODEL CALCULATION 
OF SCATTERING BY PAIRING DISORDER 

The potential which one adds to the Hamiltonian is 
then (see Shnirman et al PRE 60, 7517 (1999)). 



^1 = E(^k + ^kOn, (36) 
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where Vk = ^A(^k, and 0k is {coskx — cosky)/2, such 
that the average order parameter is Ag^k- The param- 
eter SA has dimensions of energy and is is roughly the 
amphtude of the off-diagonal modulation. The disorder 
averaged self-energies then become (with label "A" for 
order parameter modulation scattering), 



the scattering rate. As before, we examine the pole of the 
Green's function to define the approximate total quasi- 
particle scattering rate for a quasiparticle on the Fermi 
surface ^k = as 

r^((/:),cj) = -ImEo((/:^,cj) - —ImSi ((/:), cj) 
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(40) 



To gain some insight into the energy and momentum de- 
pendence of these quantities, we consider a model with a 
circular Fermi surface and approximate the d-wave order 
parameter Ak — A^ = Aq cos 20. This approximation 
will not affect the qualitative low-energy dependence of 



and estimate the <C Aq self-energies at the node (k = 
kiv, = 7r/4) and antinode (k = k^, = 0): 



riiv (^) 



4 A3 
4 Ao' 



(41) 
(42) 



and we have defined = 7rniA^o(^^A)^. It is clear that, 
as intuitively expected, the scattering of quasiparticles by 
order parameter modulations is largest near the antin- 
odes, where the mean order parameter is largest. It is 
noteworthy that this source of scattering is the only one 
considered here which gives an anisotropic in k contribu- 
tion on the Fermi surface as a; ^ which is largest at 
the antinode. 




FIG. 15: Off-diagonal scattering rate at node (squares) and 
antinode (circles) as function of energy uo for toy model where 
order parameter modulation is confined to bonds neighboring 
impurity site. 



